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ABSTRACT 


The  results  of  tests  made  with  a  two  dimensional  diffusor 
operating  with  transverse  grooves  cut  into  the  diverging  side  walls 
were  compared  with  the  results  of  similar  tests  made  using  a  smooth 
wall  diffusor.  The  results  of  this  comparison  revealed  that; 

1.  The  efficiency  of  the  ribbed  wall  configuration  exceeded  that 
of  the  smooth  wall  diffusor  for  total  opening  angles  of  approximately 
28  to  33  degrees. 

2.  The  ribbed  wall  configuration  had  an  efficiency  approximately 
120  per  cent  of  that  of  the  smooth  wall  diffusor  when  both  were  oper- 
ated at  a  total  opening  angle  of  32  degrees. 

3.  The  efficiencies  of  the  ribbed  wall  diffusor  operating  at 
large  opening  angles  varied  as  a  function  of  the  number  of  grooves  in 
the  side  wa 1  Is. 

ko      The  efficiency  increased  as  the  number  of  grooves  increased, 
but  an  unexplained  drop  in  efficiency  occurred  at  openvrrg  angles  of 
30  and  32  degrees  with  an  eight  groove  configuration. 

5.  The  maximum  number  of  grooves  tested  was  ten.  There  is  no 
evidence  to  indicate  that  this  is  an  optimum  configuration. 
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CHAPTER    I 


INTRODUCTION 


Much  work  has  been  done  and  many  papers  written  concerning  studies 
of  diffusor  performance.   Most  of  these  studies  are    limited  to  the 
study  of  diffusors  operating  with  a  total  opening  angle  of  less  than 
1£  degrees.   Cocanower,  Kline,  and  Johnston,  Ref.  \  have  proposed  a 
general  theory  for  predicting  diffusor  performance  for  a  variety  of 
two  and  three  dimensional  diffusors.  Their  theory  has  been  shown  to 
be  in  good  agreement  with  experimental  tests. 

The  studies  of  diffusor  performance  for  diffusors  operating  at 
large  opening  angles  have  not  progressed  as  rapidly  as  the  small  angle 
diffusor  work  because  of  the  complications  associated  with  the  severe 
turbulence  resulting  when  flow  separation  occurs  at  the  diffusor  walls. 
This  separation  of  flow  is  responsible  for  the  large  losses  and  poor 
performance  normally  associated  with  wide  angle  diffusors. 

Two  studies  which  deal  with  the  problem  of  wide  angle  diffusors 
are  the  papers  of  Waitman,  Reneau,  and  Kline,  Ref.  2,  and  of  Moore, 
Carl,  and  Kline,  Ref.  3.   Each  of  these  studies  deals  with  the  condi- 
tions affecting  flow  separation  in  wide  angle  diffusors. 

The  flow  separation  phenomena  observed  and  discussed  by  Moore, 
Carl,  and  Kline  in  connection  with  water  table  tests  were  also  observed 
in  the  air  flow  tests  of  this  investigation.  These  phenomena  are   dis- 
cussed in  Section  3.6  of  Chapter  III. 

Migay,  Ref.  k,  noted  an  improvement  in  the  efficiencies  of  dif- 
fusors operating  at  wide  opening  angles  when  transverse  grooves  were 
cut  into  the  diverging  walls.   His  explanation  for  this  improvement  in 


efficiency  was  that  the  "large  separation"  normally  occurring  without 
the  grooves  is  replaced,  because  of  their  effect,  by  a  series  of 
"smaller  separations."  The  resulting  losses  associated  with  these 
"smaller  separations"  are   apparently  less  than  the  loss  associated  with 
the  "large  separation."  The  efficiency  is  thereby  improved. 

Migay  cites  the  results  of  many  tests  which  support  his  claim  of 
an  improved  efficiency  when  operating  with  a  grooved  configuration. 
Unfortunately,  many  of  the  results  cited  are   not  adequately  documented. 
The  results  of  the  two  dimensional  experimentation,  for  example,  are 
presented  as  the  ratio  of  ribbed  wall  diffusor  efficiency  to  that  of  a 
smooth  wall  diffusor  operating  under  the  same  conditions.   No  informa- 
tion is  available  concerning  the  specific  conditions  under  which  the 
tests  were  conducted,  and  the  actual  efficiencies  of  each  of  the  dif- 
fusors  are    not  presented  in  the  paper. 

The  over-all  objective  of  this  study  was  to  determine  experimen- 
tally the  performance  of  a  two  dimensional  ribbed  wall  diffusor.   It 
was  hoped  that  the  optimum  operating  range  of  such  a  diffusor  could  be 
found  and  compared  with  a  smooth  walled  diffusor  operating  under  simi- 
lar conditions. 

Figure  1  illustrates  the  type  of  diffusor  used  for  the  experiment. 
The  detailed  insert  in  the  figure  shows  what  is  meant  by  the  term 
"ribbed  wall." 

Because  of  time  limitations,  the  scope  of  this  study  was  narrowed 
to  the  study  of  one  particular  rib  configuration.  The  geometry  of 
this  configuration  is  discussed  in  Section  2.h   of  Chapter  II.  A  limited 
investigation  was  also  made  to  determine  the  effect  that  varying  the 
number  of  grooves  had  upon  the  performance  of  the  ribbed  wall  diffusor. 
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ITEM  A 


TWO  DIMENSIONAL  DIFFUSOR 


DETAIL  OF  RIBBED  WALL  CONFIGURATION 


FIGURE  1 

TWO  DIMENSIONAL  DIFFUSOR  AND 
DETAIL  OF  A  RIBBED  WALL 
CONFIGURATION 
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Migay  cites  a  limited  amount  of  data  related  to  a  two  dimensional 
ribbed  wall  diffusor.  The  data  were  obtained  using  a  one  sided  dif- 
fusor  with  an  asymmetric  opening  angle  of  20  degrees.   Figure  2  shows 
the  type  of  diffusor  used.  Table  I  is  a  summary  of  the  information 
presented  in  Ref.  k   related  to  two  dimensional  diffusors. 

The  initial  rib  configuration  utilized  for  this  study  was  chosen 
to  be  approximately  equal  to  that  used  in  Ref.  1.   In  this  study,  how- 
ever, a  symmetric  diffusor  was  used. 
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TABLE  I 
SUMMARY  OF  TWO  DIMENSIONAL  DIFFUSOR  STUDIES  FROM  REF.  1 


Effect    of 

Parameter   A      (B    =  5mm 

) 

1.0 

0.857 

0.915 

1.25 

1.46 

1.64 

1.68 

1.71 

A    (mm) 

0.0 

0.500 

1.000 

2.00 

4.00 

8.00 

10.0 

15.0 

A/B 

0.0 

0.100 

0.250 

0.400 

0.800 

1.600 

2.00 

3.00 

- i 

Effect  of  Parameter  B   (A  =  20mm) 


1.61 

1.7 

1.71 

0.5 

B 

1.0 

3.0 

5.0 

unstab. 

16 

A/B 

20 

6.7 

4.0 

1.25 

Effect 

of   Nu 

mber    of 

Grooves 

N      (B   = 

:    1mm, 

A   =   20mm) 

1.04 

1.13 

1.30 

1.515 

1.63 

1.68 

1.71 

1.71 

1.71 

N 

1 

2 

3 

4 

5 

6 

7 

10 

15 

Effect  of  Number  of  Grooves  N   (B  =  3mm,  A  =  20mm) 


1.06 

1.14 

1.29 

1.40 

1.54 

1.57 

1.63 

1.64 

1.64 

N 

1 

3 

5 

7 

9 

11 

13 

15 

40 

Effect  of  Number  of  Grooves  N   (B  =  5mm,  A  =  20mm) 


0.6 

0.7 

0.78 

1.5 

1.68 

1.71 

1.71 

N 

1 

2 

3 

4 

5 

6 

15 
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CHAPTER  II 


EXPERIMENTAL  PROCEDURES 


2.1  Experimental  set-up.   The  equipment  used  for  this  study  was 
located  in  the  "New  Addition"  of  Building  230  at  the  Jet  Propulsion 
Laboratories,  Naval  Postgraduate  School,  Monterey,  California.   Figure 
3  is  a  schematic  of  the  installation.  The  source  of  air  for  the  tests 
was  the  Carrier  Air  Compressor  installed  in  Building  230. 

2.2  Pressure  readings.  All  pressure  measurements  in  the  flow 
field  were  made  using  a  Fused  Quartz  Precision  Pressure  Gaage  manufac- 
tured by  Texas  Instruments  Inc.  of  Houston,  Texas.   Prior  to,  and  twice 
during  the  time  period  of  the  tests,  the  pressure  instrument  was  cali- 
brated using  a  water  manometer.  The  calibration  remained  invariant 
throughout  the  duration  of  the  test;  and  the  value  of  0.8  inches  of 
water  pressure  per  increment,  or  count,  on  the  instrument  was  used  as 
the  conversion  coefficient  for  all  pressure  measurements.  The  details 
of  the  conversion  of  the  pressure  and  temperatjre  readings  are  dis- 
cussed in  Chapter  III,  Sections  3.1  and  3.2. 

Atmospheric  pressures  were  read  from  a  standard  wall  barometer. 
Two  pressure  readings  were  required  at  the  flow  nozzle.  The  static 
pressure  upstream,  P  ,  and  the  drop  in  static  pressure  across  the 
nozzle, Ar  ,  were  measured  and  recorded  for  the  purpose  of  determining 
the  flow  rate  through  the  system.  The  computations  involved  in  the 
determination  of  the  flow  rate  are  discussed  in  Chapter  III,  Section 
3.3. 

Two  pressures  were  measured  at  the  inlet  plane  of  the  diffusor, 
the  static  pressure  and  the  dynamic  head,  Q  .  The  static  pressure 
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was  assumed  to  be  uniform  in  the  plane  of  the  throat.  Two  static  taps 
were  located  at  the  top  of  the  diffusor  inlet  and  two  at  the  bottom. 
All  four  of  these  static  pressure  taps  were  connected  to  the  pressure 
measuring  instrument  through  a  common  lead.  The  total  pressure  was 
measured  at  three  positions  in  the  central  region  of  the  flow,  free 
of  the  boundary  layer  effect.  Figure  k   shows  the  location  of  the  pres- 
sure orifices  and  the  details  of  the  total  pressure  probe  construction, 
The  probe  was  manufactured  by  Mr.  M.  0  dell,  machinist,  Aeronautics 
Department,  NPGS. 

Figure  5  is  a  schematic  of  the  pressure  instrument  hook-up.   By 
opening  and  closing  various  pairs  of  the  six  switches,  it  was  possible 
to  measure  all  of  the  required  pressures. 

Section  3.4  of  Chapter  III  illustrates  how  the  conditions  at  the 
throat  were  calculated  using  these  pressure  measurements. 

2.3  Temperature  measurements.  Two  total  temperature  readings 
were  taken  in  the  flow  field.  The  first  was  taken  just  ahead  of  the 
flow  nozzle  and  was  used  in  the  flow  rate  calculation.   The  second 
total  temperature  measurement  was  taken  just  ahead  of  the  diffusor 
inlet.  This  temperature  was  used  in  the  determination  of  the  flow 
conditions  at  the  inlet. 

All  flow  field  temperatures  were  measured  using  a  thermocouple. 
Section  3.1  of  Chapter  III  discusses  the  details  of  the  temperature 
reading  conversions.  The  ambient  air  temperatures  were  measured  with 
ordinary  mercury  bulb  thermometers. 

2.4  Measuri  ng  and  recordi  ng  the  data .  The  rate  of  flow  through 
the  system  was  regulated  by  bleeding  the  compressor  output  to  the 
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FIGURE  k 

SCHEMATIC  OF  THE  DIFFUSOR  INLET  PLANE  ILLUSTRATING  THE  LOCATION 
OF  THE  TOTAL  AND  STATIC  PRESSURE  TAPS 
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FIGURE  5 
SCHEMATIC  OF  THE  PRESSURE  INSTRUMENTATION  INSTALLATION 
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atmosphere.  The  maximum  obtainable  flow  rate  was  approximately  2.k 
lb  /sec.  At  this  flow  rate,  the  velocity  in  the  center  of  the  inlet 
flow  channel  was  approximately  210  ft/sec. 

In  the  initial  phase  of  the  experiment  tests  were  made  at  diffusor 
angles  varying  f  rom /3  =  10°  to  /3  =  60°e  At  each  angle,  the  flow  rate 
was  varied  between  the  maximum  of  2.k   and  a  minimum  of  0.8  Ib^sec 

Tests  were  made  at  each  of  the  angles  with  smooth  walls  and  with 
ribbed  walls  configured  as  indicated  below; 

N:   number  of  grooves        =  10 

b:   width  of  grooves         =  0.125  in. 

t:   spacing  between  grooves   =  0.060  in. 

a:   depth  of  grooves         =   0.700  in. 
This  rib  configuration  was  chosen  as  a  close  approximation  to  the  rib 
sizes  used  by  Migay  in  the  two  dimensional  diffusor  study  cited  earlier 
and  summarized  in  Table  I. 

The  results  of  this  phase  of  the  experiment  indicated  that  the 
ribbed  walls  performed  best,  relative  to  the  smooth  walls,  at  /o    -   32° 
and  a  velocity  of  approximately  155  ft/sec.   Consequently,  further 
tests  were  conducted  at  this  operating  point  in  an  attempt  to  determine 
how  the  performance  varied  as  a  function  of  the  number  of  grooves.   In 
this  ph.se  of  the  experiment  the  number  of  grooves  was  increased  by 
one  for  each  successive  run.  The  results  of  these  tests  are   discussed 
in  Chapter  IV. 

The  program  used  to  perform  the  required  calculations  is  discussed 
in  detail  in  Appendix  A.   A  sample  data  sheet,  a  list  of  the  variable 
names  used  on  the  sample  data  sheet,  and  a  listing  of  input  variables 
are  included  with  the  discussion  of  the  program. 
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CHAPTER  III 


COMPUTATIONS  AND  ANALYSIS 


3.1  Temperature  measurement  conversions.  The  thermocouple 
readings,  recorded  in  millivolts,  were  converted  to  degrees  Fahrenheit 
by  the  following  equation: 

°F=  3  2.*  3  5".  98  (r^) -.455"  (mv")1- 

(3-1) 

This  equation  is  a  least  squares  fit  of  the  values  for  conversion 
from  millivolts  to  degrees  Fahrenheit  as  specified  in  the  National 
Bureau  of  Standards  Circular  5&1  for  iron  constantan  thermocouples. 
Equation  (3-1)  is  valid  only  if  the  reference  junction  is  maintained 
at  a  temperature  of  32  degrees  Fahrenheit. 

3.2  Pressure  measurement  conversions .   The  barometric  pressure, 
measured  and  recorded  in  inches  of  mercury,  was  converted  to  pounds  per 
square  inch  by  using  the  following  equation: 


^  _     ,0OOib34  t 
P-^PU-|,.OOOI8,8* 


(3-2) 

The  two  coefficients  multiplying  t,  the  temperature  in  degrees  centi- 
grade, are  thermal  expansion  coefficients  for  the  mercury  in  the  baro- 
meter and  the  brass  upon  which  the  scale  of  the  barometer  is  engraved. 
The  factor  .^912  is  the  conversion  factor  for  converting  pressure 
measured  in  units  of  inches  of  mercury  to  units  of  lbs.  per  sq.  in. 
Equation  (3-1  )>  the  conversion  factor,  and  the  correction  coefficients 
were  obtained  from  Ref.  5,  pages  1-157. 

3.3   F low  rate  measurements .   References  6  and  7  provide  general 
methods,  procedures,  and  calculations  concerning  flow  rate  measurements 
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using  flow  metering  orifices  and  nozzles.  The  flow  metering  nozzle 
used  to  determine  the  flow  rate  for  these  tests  was  constructed  in 
accordance  with  the  specifications  of  Ref.  8.  A  nozzle  was  selected, 
rather  than  a  square-edged  orifice,  because  the  calculations  were  felt 
to  be  more  easily  adaptable  to  computerized  methods. 

From  Ref.  8,  page  5,    Eq.  (la),  slightly  modified,  the  flow  rate 
through  the  flow  metering  nozzle  in  lbm/sec  is: 


w=*£A4-fti  ^|2.gc  S.CP.-PdV 


(3-3) 


The  quantities  in  Eq.  (3-3)  are   identified  below: 

W   flow  rate  (1bm/sec) 

°<   discharge  coefficient 

C*   density  coefficient 

A^  area    of  nozzle  opening  (in.  ) 

f^  area   multiplier 

\f        Reynolds  number  correction 

J^   conversion  factor  (32.17  ft.  lbm/sec2lbf) 
(tV^c\)  pressure  drop  across  nozzle  (psi) 

5^   upstream  "density"  (non-standard  units) 
By  assuming  the  flow  to  be  a  perfect  gas,  £   can  be  expressed  as 
a  function  of  the  temperature,  Tu. 


RcT, 


Re 


(3-M 
^   gas  constant  for  air  (53.35  T^   J^l  ) 
lu\   temperature  upstream  (°R) 
The  "equivalent  flow  rate,"  denoted  as  \v  ,  and  defined  by  Eq. 
(3-5),  provided  a  computational  means  of  determining  the  flow  rate. 


2k 


w  =»  — p— 


R 


9l 


H<Jc 


(3-5) 


In  Eq.  (3-5), W  has  the  units  of  in.  . 
After  substituting  for  W  and  ^   Eq.  (3-5)  becomes: 


(3-6) 


The  table  below  lists  the  values  of  the  coefficients  in  Eq.  (3-6) 
which  are  applicable  to  the  particular  nozzle  installed  for  these  tests. 

TABLE  II 
SUMMARY  OF  FLOW  RATE  CALCULATION  COEFFICIENTS 


Coefficient   Symbol 

Coefficient    Value 

Source   of   Value 

<X 

1.1 114 

Ref.   4 

£ 

Eq.     (3-7) 

Ref.    2 

ft 

Eq.    (3-8) 

Ref.    2 

<7- 

1.0 

Ref.    k 

Ai 

7.1728 

1Tr* 

If  Xs  *q —  *  then  the  coefficient  £   can  be  expressed  by  the 
equation, 


1 


5M33 


(r« 


■AZ8 


-O-x) 


^8fe 


l-.3050-^ 


\A18 


(3-7) 


If  "t^  is  the  temperature  of  the  flow  upstream  of  the  nozzle,  ex- 
pressed in  degrees  Fahrenheit,  then  the  coefficient  -f^  is  given  by  the 
equation, 


-ft  -  \+.ooi5  (~To^~) 
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Equation  (3-6)  may  be  greatly  simplified  by  combining  all  of  the 
known  constant  coefficients  into  one  constant,  c. 


C^  «<4>  Ad  fT  -  11- 30  in* 


Thus,  Eq.  (3-6)  becomes: 


\a/=  11.30  £-*V 


(P^-PaV 


P 


(3-9) 


Finally,  the  mass  flow  rate  can  be  calculated  by: 


ft* 


k 

R* 


(3-10) 


3.^  Flow  conditions  at  the  diffusor  inlet.   Velocity  and  total 
pressure  measurements  across  the  inlet  of  the  diffusor  indicated  that 
the  static  pressure  was  uniform  throughout  the  inlet  plane.   It  was 
assumed  that  the  total  pressure  at  the  inlet  was  constant  in  the  region 
outside  of  the  boundary  layer.   Based  upon  this  assumption,  the  Mach 
number  in  the  central  region  of  the  flow  at  the  inlet  was  determined 
by  solving  the  following  isentropic  relationship  for  M  : 


p,    v 


i  +  *jl  n; 


) 


Y-\ 


(3-11) 


The 


measured  values  of  p  and  p   were,  used  to  determine  M 


(3-12) 


The  static  temperature  at  the  inlet  was  determined  by  using  the 
previously  determined  Mach  number  and  the  following  equation: 


T, 


"0 


*  *j|-  <) 


(3-13) 
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Solving  Eq.  (3-13)  for  T»   yields  the  required  relationship: 

T  -       nrt\ 


\ 


'  *  V  ^  (3-1^) 

Once  1 1  ,  the  static  temperature,  was  determined,  it  was  possible 
to  evaluate  a,  the  speed  of  sound,  and  subsequently  the  midstream  velo- 
city at  the  inlet. 


(3-15) 
The  velocity  in  the  central  region  of  the  flow  is  indicated  by  ^max» 

By  assuming  the  flow  to  be  a  perfect  gas,  enough  information  was 
available  to  calculate  the  density  of  the  flow  at  the  inlet.   It  was 
assumed  that  the  density  remained  constant  across  the  inlet  plane. 

?,  *  144 


e,  -- 


R^  ^T,  (3-16) 

Equation  (3-16)  expresses  the  density  in  units  of  slugs/ft-^. 

The  Reynolds  number  based  upon  the  equivalent  diameter  of  the 
inlet  was  determined  by  expressing  the  kinematic  viscosity  of  air  as 
a  function  of  temperature.  The  expression  represents  a  fit  of  the 
data  tabulated  on  page  8  of  Ref.  9»  The  following  expression  is  con 
sidered  valid  only  for  the  temperature  range  of  50°F.  to  10A  F: 

(3-17) 
The  units  of  x)      are   ft  /sec. 

3.5   Integrated  f low  rate  at  the  dif fusor  i  nlet .   All  of  the  ca 1  - 
culations  related  to  the  conditions  at  the  inlet  and  the  performance 
of  the  diffusor  were  dependent  upon  the  total  pressure,  static  pressure, 
and  total  temperature  measurements  at  the  diffusor  inlet.   It  was  felt 
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that  some  method  of  evaluating  the  reliability  of  these  measurements 
should  be  incorporated  into  the  computational  procedures.  The  method 
decided  upon  was  to  compare  the  flow  rate  at  the  inlet,  computed  from 
the  velocity  and  density  as  determined  from  these  measurements,  against 
the  flow  rate  measured  at  the  flow  nozzle.   Since  the  two  flow  rates 
were  determined  from  completely  independent  measurements,  it  was  felt 
that  any  erroneous  readings,  or  errors  in  recording  the  data,  would 
become  apparent  when  the  two  flow  rates  were  compared. 

Prior  to  performing  the  calculations  to  determine  the  flow  rate 
through  the  diffusor,  tests  were  made  with  a  hot  wire  anemometer  to 
determine  the  velocity  profiles  at  the  inlet.  Figure  6  is  the  velocity 
ratio,  V/Vmax,  plotted  as  a  function  of  the  distance,  h,  from  the  sur- 
face of  the  wall.   Vmax  is  the  velocity  in  the  central  region  of  the 
flow.  Figure  7  represents  the  same  velocity  ratio  as  measured  from  the 
side  wall.  The  profiles  represented  in  Figs.  6  and  7  were  measured  at 
the  mid-span  position  of  the  top  and  side  walls  respectively. 

In  each  instance,  it  was  possible  to  find  analytic  expressions 
which  closely  approximated  the  velocity  profiles.  The  data  plotted 
were  taken  at  several  velocities  representative  of  all  anticipated 
operating  conditions.  Although  it  was  recognized  that  the  velocity 
profiles  vary  as  a  function  of  velocity,  tests  indicated  that  the  vari- 
ation was  minor  and  under  the  anticipated  velocity  operating  range 
would  not  exceed  the  variations  caused  solely  by  instrumentation  error. 
The  total  change  in  the  flow  rate  determination  over  the  extremes  in 
the  velocity  profile  predictions  was  found  to  be  a  fraction  0f  1  per 
cent  of  the  total  flow  rate. 
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The  expression  used  to  approximate  the  velocity  along  the  top  wa 1 

(3-18) 
When  solved  for  V  as  a  function  of  h,  the  distance  from  the  wall, 
Eq.  (3-18)  becomes: 

(3-19) 
A  similar  expression  was  obtained  for  the  shape  of  the  velocity 
profile  at  the  side  wall. 


h  -  0,60  Z 

Solving  for  V  as  a  function  of  h  yields 


-1(30(1-  v/v^" 


(3-20) 


v-v^{\-U^°-?r)' 


(3-21) 
Figure  8  is  a  cross-sectional  view  of  the  inlet  plane.   The  region 

labeled  "C"  is  the  area  outside  of  the  effect  of  the  boundary  layer. 

It  was  assumed  that  the  velocity  in  this  region,  previously  referred  to 

as  V   ,  was  constant  throughout.   Measurements  indicated  that  this 

assumption  was  reasonable. 

A  turbulence  survey  was  made  in  the  central  region  by  feeding  the 

hot  wire  anemometer  output  to  an  oscilloscope.  The  results  of  this 

study  indicated  that  the  1  maximum  velocity  perturbation  was  less  than 

1  per  cent  of  the  mean  velocity,  Vmax. 

The  regions  labeled  "F"  and  "D"  in  Fig.  8  represent  those  regions 

within  the  boundary  layer  along  the  top,  bottom,  and  side  walls 

respectively. 
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Although  it  was  not  possible  to  accurately  measure  the  velocity 
profiles  on  the  bottom  and  opposite  side  wall,  those  measurements  taken 
indicated  that  the  variations  from  wall  to  wall  and  top  to  bottom  were 
slight.  Thus,  it  was  assumed  that  the  velocity  profiles  on  opposite 
walls  were  identical. 

The  corners,  indicated  as  "E"  in  Fig.  8,  were  assumed  to  have  a 
velocity  distribution  equal  to  the  product  of  the  V/Vmax  ratios  from 
each  wall  as  a  function  of  the  respective  distances.  Figure  9  shows 
how  this  was  done.  The  velocity  assumed  to  exist  at  the  point  p, 
indicated  ir,  Fig.  9,   was  assumed  to  be  the  product  of  the  velocity 
ratio  V/Vm  x,  given  by  Eq.  (3-19)  for  the  distance  h^,  from  the  top 
wall  and  the  velocity  ratio  V/VmaX,  given  by  Eq.  (3-20),  for  the 
specified  distance,  )r\2,    from  the  side  wall. 

The  diagonal  line  indicated  as  "corner  survey  path  line"  in 

0 

Fig.  9  is  the  line  upon  which  a  survey  of  velocities  was  made.  Figure 
10  shows  the  comparison  between  what  the  above  assumption  predicted  for 
the  velocities  along  the  line  and  the  actual  measured  values  of  V/Vmax. 
The  analytic  expression  for  the  velocity  at  any  point  in  the  corner 
region  is  given  as  a  function  of  the  two  distances  hi  and  h?. 

The  purpose  of  this  investigation  into  the  velocity  profile  at  the 
inlet  was  to  obtain  a  measure  of  the  reliability,  or  validity,  of  the 
pressure  readings.  This  was  accomplished  by  using  the  measured  values 
of  the  velocity  and  density  together  with  the  estimates  of  the  velocity 
profile  to  integrate  across  the  inlet  and  determine  the  flow  rate. 
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TOP   OF    DIFFUSOR 


CORNER    SURVEY    PATH    LINE 


LIMIT   OF    BOUNDARY    LAYER 


FIGURE    9 
THE    CORNER    REGION    (E)    OF   THE    INLET    PLANE 
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The  calculation  for  the  mass  flow  rate  is  expressed  by  the  follow- 
ing equation  when  the  flow  is  considered  to  be  incompressible: 

M 


o 


(3-23) 


Table  III  is  a  breakdown  of  the  mass  flow  rate  calculations  for 
each  of  the  regions  specified  in  Fig.  9. 


TABLE  III 
SUMMARY  OF  CONTINUITY  CALCULATIONS 


Region 

Area 

Velocity 

Mass   Flow   (lb   /sec) 

Eq. 

C 

0.008 

Vmax 

(I-tt) 

D 

0.013 

Eq.    (3-19) 

(3-25) 

E 

0.020 

Eq.    (3-22) 

o,o*8  egc$  J  VWAA 

(3-26) 

F 

0.080 

Eq.    (3-21) 

°«°^f^N\^^ 

(3-27) 

The  flow  rate  determined  by  this  integration  was  compared  to  the 
flow  rate  measured  at  the  flow  nozzle  as  discussed  in  Section  3.2. 
In  general,  the  agreement  was  within  1  per  cent. 

3.6  Diffusor  performance.  The  conditions  which  existed  at  the 
diffusor  inlet  were  found  to  be  generally  steady,  and  the  measurements 
were  felt  to  be  reliable  indicators  of  the  flow  at  that  plane.   It  was 
found,  however,  that  conditions  at  the  exit  plane  were,  under  most 
operating  conditions,  extremely  turbulent.  For  large  opening  angles, 
of  particular  interest  in  this  s^udy,  the  flow  actually  reversed  direc- 
tion along  the  diffusor  walls  after  separation  had  occurred.   Moreover, 
at  opening  angles  greater  than  26  degrees,  the  flow  was  unsteady  in  that 

36 


it  would  separate  first  on  one  side  of  the  diffusor  and  then  on  the 
other.   It  was  noted  that  under  certain  conditions  the  flow  could  be 
directed  along  one  wall  of  the  diffusor  while  complete  separation  ex- 
isted along  the  opposite  wall.  The  flow  could  then  be  shifted  to  the 
opposite  wall  causing  the  separation  to  shift  sides.   It  was  also  noted 
that  the  diffusor  could  be  made  to  separate  at  both  sides  simultaneously 
causing  a  "free  jet"   flow  straight  out  of  the  inlet  throat. 

Because  of  these  non-uniform  conditions  which  existed  at  the  exit, 
it  was  not  possible  to  measure  a  single  meaningful  exit  velocity.  The 
only  flow  condition  which  could  be  accurately  measured  at  the  exit  was 
the  static  pressure,  ?2»      Because  the  diffusor  exhausted  to  atmospheric 
conditions  at  subsonic  velocities,  the  exit  static  pressure  was  taken 
to  be  equal  to  atmospheric  pressure. 

Not  being  able  to  measure  an  exit  velocity  posed  a  serious  problem 
in  attempting  to  measure  diffusor  performance.   Most  definitions  of 
diffusor  efficiency,  for  example,  require  a  knowledge  of  the  exit  velo- 
city, \l2»      Mi  gay,  Ref.  1,  defines  diffusor  efficiency  as: 

'■*       "a.  ?  (v>Vil)  (3-28) 

He  makes  no  attempt  to  explain  the  meaning  of  V2  or  ^ow  ^  was  obtained. 

A  more  convenient  and  perhaps  realistic  criteria  for  diffusor  per- 
formance are  the  coefficient  of  pressure  recovery,  defined  as 

H*        %J$t.  (3_29) 

This  measure  of  diffusor  performance  has  the  distinct  advantage 
that  all  of  the  quantities  can  be  easily  determined.   ^max»  tr,e  ve^°" 
city  of  the  flow  in  "Region  C,"  Fig.  8,  was  chosen  as  the  reference 
velocity  V ] ,  to  be  used  in  Eq.  (3-29). 
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In  order  to  compare  the  results  of  this  test  with  those  published 
in  Ref.  1,  some  scheme  had  to  be  decided  upon  to  produce  a  quantity  V2 
for  Eq.  (3-28).  One  method  used  was  simply  to  assume  that  the  central 
flow  at  the  inlet,  with  velocity  Vmax,  was  expanded  completely  and  in- 
compressibly  from  area  A]  to  area  At.  ^2  ]° s  an  ideal  expansion  of  the 
midstream  velocity  at  the  inlet. 

A  second  means  of  evaluating  the  efficiency  was  to  determine  the 
ratio  of  the  measured  to  the  ideal  coefficient  of  pressure  recovery. 

The  ideal  coefficient  of  pressure  recovery  can  be  expressed  as  a 
function  of  the  inlet  and  exit  areas, 


^  ~-   '  "fttf 


(3-30) 

Equation  (3-30)  was  derived  by  substituting  Bernoui Ill's  equation 
and  the  equation  of  continuity  for  incompressible  flow  into  Eq.  (3-29) 

3.7  Measurement  rel iabi 1 ity  indi  cat or ,   In  order  to  obtain  an- 
other indication  of  the  reliability  of  the  pressure  and  temperature 
measurements  at  the  diffusor  inlet,  a  comparison  was  made  of  the  mea- 
sured flow  rate  based  upon  the  flow  metering  nozzle  and  the  flow  rate 
calculated  by  integrating  the  velocity  over  the  area  of  the  diffusor 
inlet  as  described  in  Section  3.k.     The  quantity  referred  to  as  the 
"data  reliability  indicator"  was  defined  by  the  equation: 

Wfafot  (3-31) 
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CHAPTER  IV 
RESULTS  AND  RECOMMENDATIONS 

4.1   Results.   This  project  was  conducted  in  three  phases.   The 
preliminary  phase  consisted  of  planning,  building  and  evaluating  the 
test  equipment.  The  second  phase  consisted  of  the  development  of 
procedures  and  data  reduction  calculations.  The  final  phase  encom- 
passed the  testing  required  to  determine  the  most  desirable  operating 
range  for  conducting  further  tests  and  an  investigation  to  determine 
the  effect  of  the  number  of  grooves  upon  the  ribbed  wall  diffusor  per- 
formance. The  results  of  each  of  these  phases  of  the  project  are  dis- 
cussed individually  in  the  following  sections. 

4.1.1   Preliminary  phase.  The  comparison  of  the  flow  rate  meas- 
ured at  the  nozzle  and  the  flow  rate  computed  at  the  diffusor  inlet 
provided  a  cross  check  of  the  pressure  readings  and  the  computed 
values  of  the  Mach  number  and  velocity.  The  data  reliability  indica- 
tor defined  by  Eq.  (3-31)  reflected  the  comparison  of  these  two  flow 
rates.  The  average  value  of  this  indicator  for  the  last  50  data 
points  recorded  during  the  project  was  0.0044.  This  indicates  that 
the  agreement  was,  on  the  average,  within  0.5  per  cent.   This  is  com- 
parable with  the  accuracy  of  the  flow  nozzle  which  is  -  0.5  per  cent. 

Section  3«6  of  chapter  III  discusses  the  two  methods  used  to 
determine  the  diffusor  efficiency.   It  can  be  shown  that  the  ratio 

r^/p  is  identical  to  the  efficiency  defined  by  Eq.  (3-28).   Both 

methods  of  determining  the  efficiency  were  included  in  the  calculations 
dealing  with  diffusor  performance  since  Eq.  (3-28)  is  based  upon  the 
calculated  value  of  the  density,©  ,  whereas  the  ratio  of  pressure 
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recovery  coefficients  does  not  depend  upon  \r     .  The  comparison  of 
these  two  efficiencies  provided,  therefore,  an  indication  of  the 
reliability  of  the  measured  value  of  the  temperature  and  pressures  at 
t  he  inlet. 

The  third  indication  of  the  reliability  of  the  test  equipment 
and  procedures  was  the  repeatability  of  the  tests.  Table  IV  shows 
the  results  of  tests  run  on  separate  days  under  similar  operating 
conditions.   It  was  felt  that  the  repeatability  indicated  in  Table  IV 
was  reasonably  good  and  that  a  realistic  performance  determination 
could  be  made  on  the  basis  of  these  tests. 

The  results  of  the  tests  seemed  to  indicate  that  the  test  equip- 
ment and  procedures  could  be  relied  upon  to  an  extent  consistent 
with  the  intended  use  of  the  calculated  results.   It  was  felt  that 
sufficient  safeguards  had  been  incorporated  into  the  computational 
system  to  provide  adequate  warning  indications  in  the  event  that 
erroneous  data  were  being  processed. 

^.1.2   Deter mi  nation  of  the  most  desirable  operat i ng  range.   To 
accomplish  the  objective  of  this  phase  of  the  study,  ribbed  wall  and 
smooth  wall  tests  were  run  over  a  wide  range  of  operating  conditions. 

Table  V  summarizes  the  results  of  this  phase  of  the  tests.   Al- 
though runs  were  made  from/j  ~  1°°  to  /)  a  60°  for  velocities 
ranging  from  6$  ft/sec  to  210  ft/sec,  only  those  values  of  ft   which 
were  of  particular  interest  to  this  study  are  presented  in  Table  IV. 

Figure  11  illustrates  the  operating  range  of  particular  interest 
for  ribbed  wall  investigations.  Figure  12  shows  how  the  ribbed  wall 
efficiency  compares  to  that  of  the  smooth  wall  diffusor  in  this  region. 
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TABLE  IV 


RESULTS  OF  REPETITIVE  TESTS  UNDER  SIMILAR 
OPERATING  CONDITIONS 


fi 

Ribbed  Wa 1 1 

Smooth  Wa 1 1 

Run  No. 

V* 

Run  No. 

% 

26° 
26° 

28° 
28° 

30° 
30° 

32° 
32° 

219 
221 

219 
221 

219 
221 

219 
221 

61.75 
62.86 

59.63 
59.68 

58.94 
59.17 

56.96 
57.36 

220 
202 

220 

220 
204 

220 
207 

65.30 
66.53 

58.58 

53.51 
53.60 

47.84 
47.17 
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TABLE  V 


COMPARISON  OF  SMOOTH  WALL  AND  RIBBED  WALL 
DIFFUSOR  EFFICIENCIES-'' 


/3 

Vs 

Vr 

%/% 

20° 
22° 
24° 
26° 
28° 
30° 
32° 
34° 
36° 
38° 

4o° 

83.06 
81.21 

65.30 
58.58 
53.51 
47.84 

43.51 
39.30 
34.55 
31.66 

78.29 
69.79 
61.75 
59.63 
58.44 
56.95 
34.55 
29.82 

27.69 

0.9640 

0.9456 
1.0179 
1.1015 
1.1906 
0.7941 
0.7588 

0.8746 

•''All  data  recorded  at  a  measured  flow  rate  of 
approximately  I.78  lbm/sec.   Ribbed  wall  con- 
figuration: 

A  =  0.7  in. 

B  =  0.125  in. 
TH  =  0.06 
N  =  10 
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TABLE  V 


COMPARISON  OF  SMOOTH  WALL  AND  RIBBED  WALL 
DIFFUSOR  EFFICIENCIES^- 


/3 

Vs 

V* 

\ft* 

20° 
22° 
24° 
26° 
28° 
30° 
32° 
34° 
36° 
38° 
40° 

83.06 
81.21 

65.30 
58.58 
53.51 
47.84 

43.51 
39.30 
34.55 
31.66 

78.29 
69.79 
61.75 
59.63 
58.44 
56.95 
34.55 
29.82 

27.69 

0.9640 

0.9456 
1.0179 
1.1015 
1.1906 
0.7941 
0.7588 

0.8746 

"All  data  recorded  at  a  measured  flow  rate  of 
approximately  I.78  1bm/sec.   Ribbed  wall  con- 
figuration: 

A  =  0.7  in. 

B  =  0.125  in. 
TH  =  0.06 
N  =  10 
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4.1.3   Effect  of  varyi  ng  the  number  of  grooves  upon  r  i  bbed  wa 1 1 
diffusor  performance.  Table  VI  is  the  result  of  the  third  and  final 
phase  of  this  investigation.   Unfortunately,  time  did  not  permit  this 
phase  of  the  testing  to  proceed  beyond  N  =  10  grooves.   One  interesting 
result  of  these  tests  is  illustrated  in  Table  VI.  At  /$   ~   30°  and 

/3    =  32°,  it  was  noted  that  the  efficiency  of  the  ribbed  wall  dif- 
fusor decreased  by  approximately  20  per  cent  as  the  number  of  grooves 
was  increased  from  7  to  8.  This  phenomenon  did  not  occur  at  fj    -   30°. 
Migay,  Ref.  1,  reports  a  similar  drop  in  efficiency  for  a  three  di- 
mensional conical  diffusor  operating  at  an  opening  angle  of  31  degrees. 
His  explanation  for  this  decrease  in  efficiency  is  that  it  is  caused 
by  "acoustic  irradiation  of  inner-fin  hollows."   No  attempt  was  made  in 
this  study  to  either  explain  or  further  investigate  the  phenomenon. 

k.2      Recommendations.   Much  experimental  work  remains  to  be  done 
in  the  investigation  of  ribbed  walled  diffusors.   It  is  hoped  that 
continued  studies  in  this  field  will  make  available  data  upon  which  a 
theoretical  analysis  of  ribbed  wall  diffusor  performance  can  be  based. 

It  is  felt  that  the  next  step  in  the  investigation  of  ribbed 
wall  diffusors  should  be  an  attempt  to  obtain  data  reflecting  the 
effect  on  performance  as  other  parameters  of  the  rib  configuration 
are   varied. 

In  conjunction  with  further  experimental  studies  into  the  effect 
of  ribs  on  the  performance  of  diffusors,  it  is  felt  that  studies  into 
the  mechanism  of  flow  past  grooves  or  ribs  should  be  thoroughly  in- 
vestigated.  Fox,  Ref.  10,  has  done  much  work  in  this  field.   It  is 
felt  that  the  results  of  his  study  would  be  of  great  value  in  attempt- 
ing to  understand  the  mechanism  of  the  ribbed  wall  diffusor  operation. 


kS 


TABLE  VI 


EFFECT  OF  VARYING  THE  NUMBER  OF  GROOVES  UPON  THE  PERFORMANCE 
OF  A  RIBBED  WALL  DIFFUSOR.* 


ft 

-  28°,      A 

=    0.7,      B    = 

--   0.125,      TH 

=  0.06 

%/fc 

0.7979 

0.9339 

0.9507 

0.9706 

1.0179 

N 

6 

7 

8 

9 

10 

/3   -  30°,      A 

=   0.7,      B   = 

=   0.125,      TH 

=   0.06 

%/% 

0.7^3 

0.9587 

0.8866 

1.0525 

1.1015 

N 

6 

7 

8 

9 

10 

fi=   32°,      A 

=   0.7,      B   = 

=   0.125,      TH 

=   0.06 

yj% 

0.58^+7 

0.9505 

0.7306 

1.1378 

1.1906 

N 

6 

7 

8 

9 

10 

-All  runs  for  this  test  made  with  measured  flow  rate  of 
approximately  1.78  1bm/sec. 
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The  possible  use  of  ribs,  or  grooves,  as  a  means  of  affecting 
boundary  layer  control  in  a  variety  of  external  as  well  as  internal 
flow  situations  seems  to  indicate  that  further  research  into  the 
behavior  of  ribs  is  warranted.   It  is  hoped  that  the  equipment,  pro- 
cedures, and  experience  assembled  in  the  course  of  this  study  will 
serve  some  useful  purpose  in  any  such  subsequent  investigations. 
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APPENDIX  A 


COMPUTER  PROGRAM 


Genera  1  statements  concerni  ng  program  ' 'Diffusor  . "  A  1 1  of  t he 
calculations  discussed  in  Chapter  III  were  programmed  for  the  IBM 
360-67  Digital  Computer  in  operation  at  the  Computer  Facility,  Naval 
Postgraduate  School,  Monterey.  The  program  was  written  in  Fortran  IV 
language.  All  of  the  calculations  as  discussed  in  Chapter  III  are 
utilized  in  the  program  named  "Diffusor."  Appendix  B  is  a  listing  of 
the  program.   Included  in  the  program  listing  are  "comment"  cards  to 
help  explain  each  step  of  the  calculation.  A  basic  understanding  of 
the  Fortran  IV  language  should  be  a  1 1  that  is  required  to  follow  the 
program  operation.   In  this  Appendix  the  emphasis  is  directed  toward 
the  input  and  output  phase  of  the  program  operation. 

Program  input.  The  special  listing  of  symbols  entitled  "Program 
Input  Variables"  included  with  this  Appendix  defines  the  Fortran 
variables  names  used  as  input  for  program  "Diffusor."  Figure  13  is  a 
sample  data  page  illustrating  the  order  and  variable  names  as  recorded 
during  the  tests.  Figure  14  shows  the  field  specification  for  each  of 
the  input  variables. 

The  first  data  card  of  each  program  run  is  a  72  space  "Hollerith" 
statement.  The  first  72  spaces  of  this  card  may  be  utilized  to  include 
supplementary  or  amplifying  statements  to  be  printed  as  part  of  the 
program  output.  Any  information  punched  on  this  card  will  be  repro- 
duced directly  beneath  the  title  line  on  the  output  page. 
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PROGRAM  INPUT  VARIABLES 

MM  Number  of  sets  of  data  to  be  read 

IRUN  Number  of  run 

MO  Month  (numeric) 

IDAY  Day  of  month  (numeric) 

PATM  Atmospheric  pressure 

PINOZ  Static  pressure  ahead  of  nozzle  (  ^.) 

DLTP  Pressure  drop  across  nozzle  (AP  ) 

TN  Total  temperature  at  nozzle 

TEMP  Temperature  (°C)  for  PATM  correction 

PI  Static  pressure  at  inlet  P^ 

QI  Dynamic  head  at  inlet  (Pt-Pj) 

TTI  Total  temperature  at  inlet  (  T^i) 

BETA  Total  diffusor  opening  angle 

TPR  Outside  ambient  air  temperature  (°F) 

AA  Depth  of  grooves 

BB  Width  of  grooves 

TH  Spacing  between  grooves 

NN  Number  of  grooves 
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72 


SAMPLE  INPUT/OUTPUT  DATA  FOR  PROGRAM  DIFFUSOR 


72  Space  Hollerith  Statement  Card 


1    2 

1 

MM 

1       3 

k       6 

7       9 

10         19 

20             29 

30            39 

ko         k9 

50           59 

207 

02 

13 

29.83 

12.35 

18.95 

1.50 

22.0 

IRUM 

MO 

IDAY 

PATM 

PINOZ 

DLTP 

TN 

TEMP 

1               10 

11            20 

21             30 

31         4o 

41           50 

-0.850 

2.155 

1.50 

32.0 

57.5 

PI 

Q1 

TT1 

BETA 

TPR 

1                10 

11             20 

21             30 

31         ko 

0.70 

0.125 

0.06 

10 

AA 

BB 

TH 

NN 

FIGURE  14 

INPUT  DATA  FIELD  AND  FORMAT  SPECIFICATIONS   FOR  PROGRAM  DIFFUSOR 
(Top  line  indicates  card  spacing) 
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The  second  data  card  of  each  computer  run  will  indicate  the 
number  of  data  sets  which  are   to  be  read.  The  name  given  to  this 
"fixed  point"  variable  is  MM. 

The  remaining  data  cards  represent  the  input  to  the  program. 
Each  data  point  requires  three  cards.  As  illustrated  in  Fig.  1*+,  the 
first  of  each  series  of  three  input  cards  contains  the  following: 
IRUN,  MO,  IDAY,  PATM,  PINOZ,  DLTP,  TN,  and  TEMP.   The  second  of  each 
series  of  three  cards  includes:   PI,  QI,  TTI,  BETA,  and  TPR .   The  last 
in  the  set  of  three  cards  used  for  each  data  point  contains  the  fol- 
lowing information  related  to  the  ribs  of  the  diffusor:  AA,  BB,  TH , 
and  NN.  The  symbols  used  to  identify  the  input  variables  are   defined 
in  this  Appendix.  For  smooth  wall  runs,  this  card  was  left  blank 
resulting  in  zeros  for  all  of  the  rib  parameters. 

Special  attention  is  directed  to  the  fact  that  the  "I"  format 
variables  are  right  justified  within  their  specified  fields  whereas 
the  "F  format  variables  can  be  punched  anywhere  within  the  specified 
field  provided  the  decimal  is  punched  with  the  number.  The  field 
specifications  and  the  format  for  each  of  the  input  variables  is 
included  in  Fig.  ]k. 

Program  output .  Figure  15  is  a  sample  of  the  output  from  program 
"Diffusor."  The  same  data  used  to  illustrate  the  input  in  Figs.  13 
and  ]k   are   printed  as  output  in  Fig.  15. 

The  second  line  of  output  is  the  Hollerith  statement  read  into 
the  program  as  the  first  data  card. 

The  grouping  labeled  DATA  CHECK  is  a  print  out  of  the  raw  data 
as  read  into  the  program.   It  was  used  to  check  the  output  for 
possible  errors  in  card  punching. 


53 


oc 

I 


I 


LW 

U 

LLoo 

U.' 

t/"C 

LO 

C  C  OCCO.  OCO^C 

U-  -Jt^_  «vrv 

c 


o  »w 

ILZU. 

U)*~v. 

«v      ."V 

CCOCL 

-Jt^«- 

f  oc 

<f  r-  — 

«tCTtf 

c  cc 

•   •   • 

.-ifM»- 

«4     c  cv<c«-<f\.r^ir.o  c 

o     u  i      f~  c-'  *j  ^  «c  * 
ro    cvllil  c  ire  oc  >j^ 

a  (OC  OC      •••••• 

cr>*crc  oj  — 'c  c  c 


n      M  ii  ii  n  ii  ii  ii  ii  ii 

C  C  «  — 

i  >*r-.-<  •              — 

c  -to*  >         —  «  —           »rv 

•  CCC  tt             LL"-LL       LLLL> 

o^  •••  uj    —  sea    Kit    » 

rvj  «-ioj~  >     llc  »o-<cjj 

>-                                                                         ii  C     ao>o»-au<r<i 

<I                                                                                                       Q.  II    II    II  U      O       O      LU            LULL 

C                                                                                                S.  LL  LU      UJUJ|-J3>CC 

LL  K  a            <»UJ2CCt>->- 

^                                                                                                      >-  <  -Jtt       _J*--JLL 

*-  —     k  a^azc-K     02a 

K  L03<1  32KLL2  —  ,-h(/:l_;*-n^ 

<  ooCoc  i/iczo  xl^u-loo: 
—  <_i  «/-^-uj— -oo  iu.i<o;ra 
_j                                                                                          C  ZU.3  LLt-U>-CG— 2      »-  — 

>                                                                      •  k_j  ac    «-c     cci/i>->- 

u.                                                                        •-<  lu2ll  zxwZI'-u.u.'lc 

i-u.  ac<z>ui-i-eitz? 

oj                                                                         11  <jC  Ctizu-ii.«<au.a 

u                                                                        2  a<u.  >oaZa— 1     •-•- 

2                                                                                                »-  >QC  •►-                         3-iOC 

<  »i—  LL  llll  <o<^-«-< 

I  G3LL  LL_J                         LL_JLLLLLL 

o                                                                            c  JCu.  cz                a«acu.u 

C  co  ••u.lloc  #.<_>•—  <o>-lulu 
a 

a 


ex 

C  a 

«/;  C 

~  LO 

li  ~ 

LL  LU 

>—  LL 

c        — 

c 

_l 

_l        I 

<        < 

1  or 
C 

c        c 

LL  (X 

cc  a 

cc 

—  a 

a  c 

LL 

—1 

<s         < 

2  h- 
(-.           < 

>—       a 

Z  K- 

LU  3 

s.        a 

•-        1- 

O  3 

O 

C  v. 

h-  3 

a 

2 


C 

-JCLL 

C2 

LO 

••u.Luar 

..<j>»- 

• 

K 

K 

<->r- 

"^ 

3 

O     •    •    • 

inir> 

a 

a 

LL222 

a  11 

K 

h- 

C«->*-*- 

•a 
ceo. 

>-ih- 

c 

c 

CJC  C  C 

L..LCOC 

11 

r-  (\jC  vT 

0 

LL 

c  —  r-v 

\- 

_j 

•  •  •  • 

-JO 

rv, 

OJC.C  c  c 

Clt 

k~ 

LL 

r«"             — 1 

• 

C 

<  < 

•— <<< 

.-j 

2 

— — c 

on  LOLL 

ai/n/jo' 

II 

2 

a  Q.C 

q  a  a 

II  II  II  II  II 

r-l 
Cl- 

C 
—J 

c 

•— ■ 

LOK 

LL 

C  f^-Ovj 

c  10  f-  e 

1- 

CO 

r-r-^ 

c  Lor^c 

~- 

• 

I 

ir -to 

>Or-rOc_ 

<\J 

O 

crcoa 

•  lO>£<£i 

0" 

■H 

Z 

•  •  • 

nI       •      •      • 

LU 

II  LO 

►— • 

■tC't 

.-l^^^ 

—        >  — 

rv,LO 

-*      OD 

t— ♦•— <p-( 

LU— —  CZ 

O-i 

OJ 

II 

_ior«c  — 

2    • 

00 

II    II    II 

II  II  II 

C a 

-<(M 

0 

z        eu~. 

a 

• 

LU 

•— — 

<</".«/••       LU 

11 

ro 

a. 

000— '— 

LULUZ  > 

vH 

2 

<J 

LL  •—<(_) 

a  >>LUO 

O 

I 

LL»- 

h- 

a  ►-»—►- 

C  C  C  LL  O 

y- 

a<r 

< 

a<ci- 

wOCJn 

f 

wtt 

Cj 

t-K<l 

ZaakC 

c 

3 

LOLL 

OLO      »- 

LLLDCLLI 

rn 

—  </lO. 

LL 

•-     l-oo 

LL            COLL 

ceo 

0 

LLUJZ 

O 

QCI— LL 

—  LLLL       C 

•  LO 

LU 

-J  a  ll 

2 

LLLU_JI- 

DOCC 

O-OT 

a 

rvja»- 

< 

I_J2  — 

w           Zet 

oj   • 

3 

KJ 

2 

a2»— x 

XI  —LL 

HC 

00 

C<I-J 

CL 

O0i—»»UJ 

•  • 

<v-  >~  uer 

z  1 

< 

2l-< 

C 

c—    ~ 

«2 

i-oa<2 

¥-  II 

LL 

~_JI- 

LL 

2.       -< 

c 

LU»— LUQ.  ^ 

<lf-i 

X 

—  LUO 

a 

K-<l-OJ 

"- 

(T3C^Z 

a.  a 

••O.OI- 

LL 

t.<aaa 

k- 

* 

» 

h- 

a. 

K 

< 

O 

LU 

3 

3 

Of 

LU 

»— 

a 

a. 

CL 

3 

I 

< 

Z 

0 

z 

c 

O 

CC 

*— 

00 

1-^ 

— 

3 

LL 

< 

2 

LL 

Z 

K 

D 

LL 

C 

< 

-1 

— 

U 

O 

LL 

O 

ro 

t. 

c 

< 

• 

1 

CC 

II 

O 

CO 

3 

** 

Ll. 

c 

Ll. 

LL 

i-< 

I- 

a 

< 

2: 

-J 

LA 

< 

^ 

1 — 

cc 

LJ 

C3 

-J 

LU 
CC 

O 
CC 

^ 

13 

D- 

l; 

C3 

LL 

i—t 
Ll. 

fe 

<I 

CL 

1- 
3 

a 

O 

3 

LU 

tr 

_l 

_i 

21 

a 

< 

^ 

00 

0 

LL 

a 

00 

< 

LL' 

2 

a 

LL 

C 

h- 

h- 

< 

<I 

O 

U 

*— 

3 

c 

C 

z 

_i 

-J      c 


< 

a 

2 

<N  13  I  <        Z  O       2  LU 

000  C  o  a     ^-  00     *-  ct 

O^' 

2 

31  C  <  _l  —  < 

at  u  o  ll  a  o 


ik 


The  input  sections  under  the  headings  of  "Flow  Rate"  and 
"Diffusor  Performance"  are  the  converted  values  of  the  input  data. 

The  output  sections  under  these  two  headings  are  those  values 
calculated  in  accordance  with  the  discussion  in  Chapter  III. 
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APPENDIX   B 


FORTRAN     IV    G    LEVFL     1,     MOD    0 


MAIN 


DATE    =    68^61 


^5/12/3' 


OC02 

COCA 

c  o  n  5 

OOC6 

CCC7 
C0C8 


OCu<3 
CCl^ 
0111 
0012 
0")13 
0014 
CO  15 
C316 
CC17 
0  C  1  8 

on" 

0320 
0*21 
0122 
00  2  3 
H024 
C>2  5 
"026 
'J  027 


002? 
0029 
t«J?c 
00  31 
1032 
00  33 
003* 

C136 
JO  37 


c  :3P 

0039 
004u 
3C41 
3042 
00  A3 
0C44 


0045 


0046 


G047 


;ar 


ASSIGN  VALUES  TO  FIXED  CONSTANTS  NEEDED  FOR  COMPUTATIONS 

TSTD=5ie.7 

PSTD=29.92*.4921 

TEXAS  =  .P 

G=  32.17 

R  =  53.35 

GAM=  l.A 

MM  =  0 

A=2.*G*778.3 

READ  IN  INPUT  DATA  ANO  PRINT  OUT  FOR  DATA  CHECK 

READ  (5,lOrP) 
RFAD  (5,lnl  NUM 
15  READ  (5,20)  I  RUN, MO, I  DA Y, PATM, P1N07 . DLTP , TN , TEMP 
READ  (5,25)  P 1 , Ql, TT 1 , BETA ,TPR 
RFAO  (^,26)  AA.RR,TH,NN 
MM=MM+1 
WRITE  (6,399) 
WRITE  (6,4^0) 
WRITE  (6,n00» 
WRITE  (6*4011  IRUN,MO,IOAY 
WRITE  (6,4n2) 
WRITE  (6,4^3)  RETA 
WRITE  (6,4~4)  RR 
WRITE  (6,4->5)AA 
WRITE  (6,4(6)  TH 
WRITE  (6,4<~'7l  NN 
WRITE  (6,380) 

WRITE  (6,381)  PATM,P1N0Z,DLTP,TN,TFMP 
WRITE  (6.382)  PI ,01 , TT1 ,TPR 

CONVFRT  DATA  INTO  D IMFNS I ONALLY  CORRFCT  UNITS  OF  mfaSMRF 

PATM=PATM*(l.-(  (.0Cni634*TEMP)/(l.  +  ."'*'MR18*TFMPl  l)*a4<>12 

t>2  =  PATM 

DLTP=DLTP*TEXAS*.03613 

D1NOZ=PATM+P1NOZ*TEXAS*.0  3613 

Q1=Q1*TE*AS*. 03613 

Pl=PATM  +  Pl*.03f«13*TEXAS 

PT1=P1*01 

TT1=  3  2.*3  5.98*TTl-.435*TTl**2+459.7 

TN=  32.*35.<58*TN-.435*TN**2 

TNR=TN*459. 7 

REGIN  FLOW  RATE  CALCULATIONS  USING  FLOW  NQZZLF  RFADIMGS 

X=DLTP/P1N07 

Y=( l.-X)**1.42P 

Z=( l.-X)**T.2Pf 

EPS=S0RT(2.433*Y/X*(l.-Z)/( l.-.30  5*Y» ) 

FT=l.+.0COp15*(TN-6r.  ) 

WE0IV=11.3*EPS*FT*S0RT(X) 

W  =  WE0IV*P1N07*S0RT(G/(R*TNR1  ) 

COMPUTE  MACH  NUMRFR 

AMACH=S0RT(2./(GAM-1.I*((PT1/P1)**( (GAM-1. ) /GAM, -I. ) ) 
CCMPUTF  STATIC  TEMPERATURE  AT  INLFT 

Tl=TTl/( l.+t GAM-1. )/2.*AMACH**2) 
CCMPUTE  VELOCITY  IN  THE  CFNTRAL  FLOW  REGION 

VMAX  =  AMACH*S0RT(GAM*R*G*T1  I 
COMPUTE  THE  DENSITY  OF  THE  FLOW  IN  THF  CENTRAL  REGTON 

R0  =  P1*144. /(R*G*T1  ) 
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FORTRAN     IV    G    LEVEL     U     MOO    n.  MAIN                                         DATF    ■    68061                           r>5/12/30 

C  RFGIN    THE    NIIMFRICAL     INTEGRATION    TO    DETERMINE    THE    FLOW    PATF    AT     THP 

C  INLET    PLANF 

C 

C049  WSUBD=0.0 

C050  WSUBE=  r  ,<-• 

0051  WSUBF=  0.0 

00*52  H  =  -.C25 

0C53  00  IOC  N=ltll6,5 

0054  h=H+.05 

00*5  5  V  =  VMAV*(  l.-.nif67*(  ALnr,(1.2/H)  »**?» 

0056  103  WSUBF=RO*G*V*.f  5*4. 8/144. +WSUBF 

C057  WF  =  2.*WSIIBF 

0'^58  WC=R0*G*VMAX*4. 8*1.6/144. 

0059  H=-.025 

0060  OH  20^    N=l,56.5 
OCfcl  F=H*.05 

CO  62  V  =  VMAX*(  1. -.03333*1  AL0GI.60/H)  )**?) 

0063  WSUB0  =  R0*G*V*.i*5*1.6/144.*WSUBD 

0064  Hl=-.025 

0065  00  2C"  J  =  1.116,5 
C!66  H1=H1*.05 

00  67  V=VMAX*(  l.-.01667*(AL0G<1.2/Hl>  I  **2  1  *  (  1  .  -  .  ^33  3  3*  (  ALDG<  ,6/H) )**?) 

0C6R  2)0  WSUBE=RO*G*V*. 0025/144. ♦  WSUBE 

006^  l*E  =  WSUBE*4. 

«C*C  wn=?,*wsuBO 

HC71  WCALC=WC+WO*WF*WF 
C 

C  COMPUTE  THE  COFFFICIFNT  OF  PRESSURE  RECOVERY 
C 

0072  CPR=<P2-P1>/01 
C 

C  CCMPUTF  THF  APFA  RATIO  OF  THE  OIFFUSOR  FOR  THE  SPECIFIED  OPENING 

C  ANGLE 

C 

0073  BETA=BETA/57.3 

0074  ALFA=RFTA/?. 

0075  BPTA=BFTA*57.3 

C076  ARATIO=. 33333*<3. 1875*16. 1*TAN( ALFA)-. 1875/COSIALFA)  > 

C 

C  COMPUTE  REYNOLDS  NUMBER  BASED  UPON  MIOSTRFAM  VFLOCITV  AND  THF 

C  HYDRAULIC  DIAMETFR  OF  THE  INLET  FLOW  CHANNFL 

C 

0?77  VIS  =  .nOO'>O"'6*(TTl-45Q.7l  +  .00012 

0"!78  RE  =  VMAX*.4/VIS 

C 

C  COMPUTE    THE     IDEAL     COEFFICIENT    OF    PRESSURF    PFCOVERY 

C 

0C7S  PRSTAR=l.-(l./ARATIO)**2 

C 

C  COMPUTE  THE  PFPFPRED  FLOW  RATE  (NORMALIZED  TO  STANDARD  DAY 

C  CONDITIONS) 

C 

0C80  TPR=TPR+459. 7 

CC81  WSTO=W*SQRT(TPP/TSTDJ/IPATM/PSTD) 

C 

C  CCMPUTF  THF  VFLOCITV  AT  THF  EXIT  BASED  UPON  AN  IOFAL  FXPANSION 

C  OF  THE  MIO  STRFAM  FLOW 

C 

0182  V2=VMAX/ARAT10 

C 

C  CCMPUTF  EFFICIFNCY  BASED  UPON  THE  DEFINITION  OF  FO(3-28) 

C 

0  0  83  EFNCY=<P2-Ol)*144./(.5*RO*<VMAX**2-V2**?) ) 

C 

C  COMPUTE  EFFICIENCY  BASED  UPON  THE  IOFAL  COEFFICIENT  OF  PRFSSIIPP 

C  RFCOVERY 

C 

0;«4  FFFN2=  CPR/PRSTAR 

C 

C  OFTERMIKiF  THF  DATA  RELIABILITY  INDICATOR 

C 

CJ85  ERR=  1.0  -(W/WCALC) 
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FORTRAN 

CC86 

0087 

ones 

0C89 

0C90 

UC<31 

CO<52 

0093 

0C94 

cnqs 

0096 

CC97 

0098 

0C99 

0100 

C1C1 

OH  2 

0K3 

01C4 

OK  5 

01C6 

0107 

01C8 

C1C9 

0110 

0111 

0112 

0113 

0114 

0115 

0116 

C117 

0118 

0119 

012C 

0  121 

0  122 

0  123 

0124 

012  5 

IV  R  LEVEL  1.  MOO  0 


MAIN 


DATE  =  68061 


"5/12/3" 


PRINT  ALL  VALUES  AS  COMPUTED 


0126 

0127 

0128 
C129 

0130 

0131 

0132 

C133 
0134 
0135 
0136 
0137 
0138 
0139 
C140 
01^1 

0142 

0143 


WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 
WRITE 


(6,4">8 
(6,409 
(6,410 
(6,411 
(6,412 
(6,413 
(6,409 
(6,414 
(6,415 
(6,416 
(6,417 
(6,418 
(6,419 
(6,420 
(6,421 
(6,423 
(6,424 
(6,425 
(6,426 
(6,427 


10 
20 
25 
26 
380 
381 


IF     (MM.LT.NU 
FORMAT     (I?) 


FORMAT  (313, 
FORMAT  (5F10 
FORMAT  (3FK 
FORMAT  (15X, 
PORMAT  (25X, 
1  TN=',F5.2, 
382  FORMAT  (25X, 

15.2//) 
399  FORMAT  (1H1/ 
40"  FORMAT  (30X, 
1ALUATI0N  OAT 
411  FORMAT  (15X, 
402  FORMAT  (15X, 
433  FORMAT  (25X, 

404  FORMAT  (25X, 

405  FORMAT  (25X, 

406  FORMAT  (25X, 

407  FORMAT  (25X, 

408  FORMAT  ( 15X, 

409  FORMAT  (  20  X  , 

410  FORMAT  (25X, 
1SS)     =',F7 

411  FORMAT  (25X, 
10W  RATE=' ,F7 

412  FORMAT  (25X, 
1  RATE    =  ',F 

413  FORMAT  (15X, 

414  FORMAT(25X,' 
1SSURE  RECOVE 

415  FORMAT  (25X, 
IIONS: •) 

416  FORMAT  (25X, 
1ENTRAL  CORE) 

417  FORMAT  (25X, 
KCENTRAL  COR 

418  FORMAT  (76Xf 

419  FORMAT  (76X, 

420  FORMAT  (71X, 

421  FORMAT  (71X, 

423  FORMAT  (71X, 

424  FORMAT  (71X, 

425  FORMAT  (71X, 

426  FORMAT  <//15 

427  FORMAT  (25X, 
1.4) 

1J0O  FORMAT  (72H 
I 
END 


) 
) 

)  P1N0Z,W 

)  DLTP,WFOIV 

)  TN,WSTD 

) 

) 

)  PATM,CPR 

)  PI 

)  PT1,VMAX 

)  P2,R0 

)  RE 

1  AMACH 

)  ARATIO 

)  WCALC 

)  PRSTAR 

)  EFNCY 

)  EFFN2 

» 

)  ERR 

M)  RO  TO  15 

5Fn.4) 

.4) 

.4,110) 

•DATA  CHECK'/) 

•PATM=',F6.3, •       P1N07=», 

•       TEMP=« ,F5.2,'   DFR  C 

•P1=',F6.3,«      Q1=',F6.3,' 


P6.3,'       OLT 
) 

TT1=« ,F^.2 


////) 

•TWO  DI 
A'//) 
4HPUN  ,13 
•CONFIRUP 
•PFTA  (DI 
•WIOTH  OF 
•DEPTH  OF 
•SPACING 
•NUMBER  0 
•FLOW  RAT 
•  INPUT: • , 
•P1(N07?L 
.4,"  LP/ 
•DELTA  PR 
.4,'  SO. 
•TOTAL  TF 
7.4, •  LR 
•DIFFUSOR 
ATMOSPHER 
RY  =  ',F9 
•PI  (INLF 


MENSIONAL  RIBBED  WALL  <M 
I2,1H-,I2,3 


,80X,5HDATE 
ATION: •/) 
FFUSOR  ANRLE) 

RROOVES  (B) 

RROOVFS  (A) 
BETWEEN  RROOVFS 
F  RROOVES  (N) 
E,  MEASURED  WITH 
30X, 'OUTPUT 
E)  =• 

SEC  ) 

ESSURE     =',FR.4, 
IN.'  ) 

MPERATURE  =',F8.4, 
/SEC'//) 

PERFORMANCE  DATA'/) 
IC  PRES.   =',F8.4, 
.4) 
T  STATIC)  =',F8.4, 


(T 


:• ) 

,F8.4, 


•PT1  (INLE 
=',F9. 
•P2  (EXIT 
E)   =',E9. 
•REYNOLDS 
•MACH  NO. 
•AFEA  RATI 
•CALCULATE 
•IDEAL  PRE 
•EFFKIENC 
•EFFICIENC 
X,  'OATA  RF 
•1.0  -  (FL 


T  TOTAL)  =' ,F8.4, 
4,'   FT/SEC) 
STATIC)   =',F8.4, 
4,'   SLURS/CU  FT' 
NO.  (EOV.OIA. )=', 
(  INLET  CORE)   =', 
0  (FXIT/INLET) 
D  MASS  FLOW  RATE 
S.  RECOVERY  COEF. 

Y  (USINR  IDEAL  V? 

Y  (PR/PR  IDEAL) 
LIARILITY  INDTCAT 
OW  RATE  MEASURED/ 


=  •  ,F 
=  ',F 
=  ',F 
)  =  ',F 
—   • 

.022  I 
PSI 
PSI 

DFO. 

PSI 

PSI 
PSI 
PSI 


FFUSOR  P 

H-68//// 

10.4,' 
10.4,' 

K.4,' 
10.4,' 
.  13//) 
NCH  FLOW 

A',13X,« 

A',13X, • 

F',l?Xf 

',11X,  • 
A', 12*,' 
A'.12X, • 
A',12X,  • 


P=',F7.3.» 

,  •      TPR=',F 

ERFORMANCF  EV 
/) 


DFR'  1 
IN.  •  ) 
IN.') 
IN.  '  ) 

NP?ZLF'  /) 

FLOW    RATF     (MA 

FOUTVALENT    FL 

•RFFEPFO    FlOW 

COEF.     OF    PRF 

INLFT  CONPIT 

VMAX  (C 

DFNSITY 


E9.4) 
FP.4) 

=',F9.4) 

=  ',F9.4,'   LB/SFC) 

=',F9.4) 
)  =',F9.4) 

=  • ,F9.4) 
OR'/) 
FLOW  RATE  CALCULATED)  =»,F7 
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